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FOREWORD 


Tree  measurement  data  are  continually  being  collected  for  research 
and  administrative  purposes.  Once  on  hand,  these  data,  to  be  useful, 
must  be  reduced  before  they  can  yield  meaningful  results.  The  compila- 
tion process  is  often  arduous  and  time  consuming.  Presented  here  are 
ways  to  improve  and  speed  up  this  task. 

High-speed  computers  are  in  common  use  for  certain  forestry  opera- 
tions. Undoubtedly,  such  computers  will  be  used  with  greater  frequency 
in  the  future..  Nevertheless  there  are  many  places  and  instances  where 
desk  calculators  are  being  used  and  will  continue  to  be  used.  Both 
methods  for  compiling  data  have  been  considered. 
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COMPILATION  OF  TREE  MEASUREMENT  DATA 
BY  HAND  OR  COMPUTER 


By 

Paul  M.  Haacki/ 


SUMMARY 

A method  for  compiling  tree  measurements  by  hand  or  computer  is 
presented  along  with  the  work  form  and  procedures  for  its  use. 

A part  of  the  form  is  used  to  tabulate  the  necessary  identification 
and  measurement  data  from  which  IBM  cards  can  be  punched  and,  in  turn, 
processed  by  a computer.  The  chief  output  of  the  FORTRAN  program,  written 
for  the  IBM  1620,  is  cubic-foot  volume  by  Smalian's  formula  for  specified 
portions  of  the  stem,  and  International  1/4-inch  and  Scribner  board-foot 
volumes  for  sawtimber-size  trees. 

The  cross-sectioned  part  of  the  form  can  be  used  for  a graphical 
solution  of  cubic-foot  volume  and  the  estimation  of  the  diameter  and  bole 
length  at  any  point  on  the  tree. 

The  computer  program  will  need  little,  if  any,  adjustment  to  be 
usable  elsewhere  when  any  of  the  three  volume  classes  is  desired. 


1/  Research  Forester  and,  since  1961,  the  biometrician  for  the 
Experiment  Station. 
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THE  PROBLEM 


Tree  volume  tables  or  equations  for  some  areas  and  species  are  either 
nonexistent  or  could  benefit  from  revision  because  of  changing  standards 
and  inadequate  samples.  Also,  particularly  when  no  equation  is  developed, 
^differences  between  diameter  limits  and  practices  of  obtaining  stem 
measurements  on  various  projects  result  in  tables  not  specifically  appli- 
cable to  each  other's  needs.  Much  past  data  are  directly  usable;  others 
require  partial  adjustment  by  arithmetic  or  graphical  means;  with  some 
there  is  uncertainty  as  to  butt  swell  or  inconsistencies  in  measurements. 
Finally,  to  facilitate  incorporating  existing  tree  measurement  data  with 
that  being  currently  obtained,  it  is  desirable  to  standardize  the  record- 
ing and  compilation  of  all  such  information. 


SCOPE 

The  Tree  Measurements  form  (fig.  1),  which  is  described  below,  is 
well  suited  for  longhand  computation  of  cubic-foot  volumes  for  specified 
portions  of  the  stem,  either  graphically  or  by  Smalian's  formula  or  any 
other  similar  formula.  Hand  methods  for  determining  board-foot  volumes 
by  16-foot  logs  are  done  by  table  look-up  or  by  formula  for  results  compa- 
rable to  computer  output. 

The  computer  program  (fig.  2)  solves  cubic-foot  volume  inside  bark, 
using  Smalian's  formula,  for  four  levels  of  utilization.  Board-foot 
volumes  are  computed  from  equations  for  the  International  1 /4-inch  rule 
and  Scribner  rule  in  16-foot  logs  to  the  closest  log.  Cubic-foot  ratios, 
board-foot/cubic-foot  ratios,  form,  and  the  number  of  sawlogs  are  also 
found.  Descriptive  information  concerning  the  tree,  too,  is  part  of  the 
output . 

Complete  generalization  of  the  computer  program  is  impossible.  Yet, 
for  example,  few,  if  any,  changes  are  needed  when  utilization  limits 
differ  from  those  herein  described.  Re-defining  the  intent  of  some  of 
the  codes  often  will  suffice  to  broaden  the  program's  utility. 


COMPILATION  METHODS 

Field  data  are  recorded  on  the  Tree  Measurements  form.  The  small 
block  at  the  top  of  the  form  contains  descriptive  information  pertaining 
to  the  tree  which,  if  machine  computation  is  employed,  will  be  duplicated 
in  the  first  35  columns  of  the  IBM  cards  for  that  particular  tree.  Also, 
if  machine  computation  is  used,  columns  1,  2,  3,  6,  and  7 in  the  large 
block  on  the  form  will  be  punched  in  card  columns  36  through  49.  Thus, 
each  row  in  the  large  block  for  which  entries  are  made  will  require  one 
IBM  card;  columns  1 through  35  of  the  IBM  card  contain  descriptive  material 
for  the  tree;  columns  36  through  49  contain  the  specific  information  for 
the  section  of  the  tree. 
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Double  bark  thickness  (coded  D.B.T.  in  column  6 of  the  large  block) 
is  not  used  in  the  FORTRAN  program.  However,  by  punching  it  and  diameter 
inside  bark  (coded  D.I.B.  in  column  7 of  the  large  block),  bark  ratios 
could  be  found  later. 

Note  that  the  decimal  places  are  implicit  and  are  not  key-punched 
in  the  cards. 


Tree  Measurement  Form  Entries  Described 

Following  is  a description  of  each  entry  on  the  Tree  Measurements 
form--its  definition,  intent,  and  IBM  coding. 

Upper  small  block. — The  portion  of  the  block  above  the  double  line 
lists  descriptive  information  for  the  tree  which  is  duplicated  on  all 
cards  for  that  tree  (one  card  for  each  tree-section) . Information  below 
the  double-line  is  not  coded  for  IBM  use. 


Nomenclature 

: IBM 

coding 

: 

on  form 

• Name  : 

Digits—^ 

• Columns 

Definition 

Tree  No. 

NOTRE 

XXX 

1-3 

Tree  number. 

Area 

Left  part  Ell 

XXX 

4-6 

Left  digit — subregion; 

Spp. 

Right  part  Ell 

XX 

7-8 

right  two  digits-- 
unit . 

Species  code.  Hard- 

Age 

Left  part  EI2 

XXX 

9-11 

woods  coded  70  or 
over;  softwoods  under 
70. 

Age  at  stump. 

Type 

Right  part  EI2 

XX 

12-13 

Local  forest  type. 

Elev. 

Left  2 digits 
of  El 3 

XX 

14-15 

Elevation  in  hundreds 

FCC . 

Next  2 digits 
of  EI3 

XX 

16-17 

of  feet. 

Form  class  code: 

01  DIB  § Ht.  above  B.H. 
DIB  at  B.H. 


02  DIB  \ Ht.  above  swell 
DIB  above  swell 


-3- 


Continued 


Nomenclature 
on  form 

FCC . 


CuC . 
Bd.  C. 


D.B.H. 

Height 

FN. 

FD. 


IBM  coding 

Name  : Digits— ^ I 

Next  2 digits 

of  EI3  xx 


[ Definition 

Columns  : 


16-17  Form  class  code: 

03  DIB  @ top  first  log 
DOB  at  BH 


04  DIB  @ top  first  log 

DOB  above  swell 

05  (Specify  numerator) 
(Specify  denominator) 


Next  2 digits 


of  EI3 

XX 

18-19 

Cu.  ft.  vol.  deter- 
mination code: 

01 

Smalian ' s 

02 

Graphic 

03 

(Specify) 

Right  2 digits 
of  EI3 

XX 

20-21 

Bd.  ft.  vol.  deter- 
mination code: 

00 

None 

01 

International  1/4- 
inch 

02 

Scribner 

03 

Both  Int.  and 
Scribner 

04 

(Specify) 

DBH 

x^c 

22-24 

Diameter  breast  high 
o.b.  (4-1/2  ft.  above 
ground  on  uphill  side) 
or,  in  case  of  swell- 

butted  trees,  D.O.B. 
1 foot  above  swell. 


HT 

FN 


FD 


xx^x 

x^c 


x: 


25-29  Total  height,  ground 

to  tip. 

30-32  Value  of  numerator  of 

form  class  ratio  given 
in  FCC. 

33-35  Value  of  denominator 

of  form  class  ratio 
given  in  FCC.  FD>0. 
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Continued 


Nomenclature  : 

IBM 

coding 

Def i ni t i nn 

on  form  : 

Name 

• Digits.!/  : 

Columns  : 

Location 

-- 

-- 

Write-in  for  locality 

Field-- 
date,  by 

-- 

-- 

Write-in  for  field 
party  and  date. 

Of f ice — 
date,  by 

-- 

-- 

Write-in  for  office 
checker . 

1/  All  decimal  points  implicit--none  punched  in  IBM  cards. 

Upper  large 

block . --This 

portion  of 

the  form  lists  specific  entries 

for  each  section  of  the  tree  from  ground  to  tip.  If  tree  volume  compi- 
lations are  to  be  done  on  the  computer,  only  columns  1,  2,  3,  6,  and  7 
are  keypunched  in  IBM  cards.  No  further  columns  need  be  considered  un- 
less hand  calculations  are  made.  Each  row  requires  an  individual  IBM 
card. 

Nomenclature  : 

IBM 

coding 

Definition 

on  form  : 

Name 

’ Digits^/  • 

Columns  : 

Section 

— 

— 

— 

-- 

No. 

NOS 

XX 

36-37 

Sequences  numbering 
of  section. 

Kind 

KINDS 

XX 

38-39 

Labels  section  as  to 
its  location  in  stem. 
For  example  

01  Interior  Alaska--sec tion  reaching  1-ft.  stump  (or  as 

specified) . 

Coastal  Alaska--section  reaching  stump  equal  in  height  to 
D.B.H.  but  not  under  2 ft.  or  over  4-1/2  ft.  (or  as 
specified) . 

02  Interior--sec tion  reaching  2-ft.  stump  (or  as  specified). 
Coastal--open. 

03  Open. 

04  Interior  and  Coastal--section  (s)  reaching  D.B.H.  or  D.O.B. 
above  swell. 
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Continued 


Nomenclature 
on  form 


05 


IBM  coding 

yi * Definition 

Name  : Digits-t'  : Columns  : 

2 / 

Interior  and  Coastal  PT— ' --sec tion  (s)  above  D.B.H,  but  not 
reaching  8.0-,  6.0-,  or  4.0-inch  D.I.B.  tops  or  the  tip, 
having  lengths  not  over  8.15  ft. 

Interior  and  Coastal  ST.^/-- sect  ion  (s)  above  D.B.H.  ter- 
minating  at  the  ST  tops— 7 in  lengths  not  over  8.15  ft. 
such  that  log-intervals  of  16.3  ft.  accumulate  from  the 
stump  (coded  KINDS  "01")  as  the  ST  top  is  approached  or 
reached.  Even  if  the  6.0-  or  8.0-inch  D.I.B.  points 
occur  at  or  below  the  ST  top,  code  the  section(s)  reach- 
ing it  KINDS  "05"  to  get  cubic-foot  volume  to  this  point. 


06 


07 


08 

09 


Interior  and  Coastal--sec tion  (s)  reaching  the  8.0-inch 
D.I.B.  point  in  PT,  or  in  ST  where  this  point  occurs 
above  the  ST  top. 


Interior  and  Coas tal--section  (s)  reaching  the  6.0-inch 
D.I.B.  point  in  PT,  or  in  ST  where  this  point  occurs 
above  the  ST  top. 


Interior  and  Coastal--section  (s)  reaching  the  4.0-inch 
D.I.B.  point  in  PT  or  in  ST. 

Interior  and  Coastal--the  section  (just  one)  reaching  the 
tip . 


Length (ft.)  SLGTH  xpcx 


40-43  Length  of  section  in 
feet  and  hundredths. 
Two  decimal  places. 


Cumu 1 . 
Length(ft.) 


D.O.B. (in. ) 


D.B.T.(in.)  DBT 


Cumulative  length  of 
sections  to  this  point. 
Total  for  tree  is  also 
entered  in  SMALL  BLOCK 
item  "Height".  Two 
decimal  places. 

Diameter  outside  bark 
at  the  top  of  the  sec- 
tion. One  decimal 
place . 

44-46  Double  bark  thickness 
at  the  top  of  the  sec- 
tion. One  decimal 
place. 
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Continued 


Nomenclature  : 

IBM  coding 

on  form  : 

Name 

• Digits.!/ 

• Columns 

1—s  v-  oL  J-  1 1 J-  L>  -L  LJ  1 1 0 

D.I.B. (in.) 

DIB 

47-49 

Diameter  inside  bark 
at  the  top  of  the  sec- 
tion. One  decimal 
place. 

(D.I.B.)2 

(sq.  in.) 

Cross-sectional  area 
in  sq.  in.  at  top  of 
section.  Two  decimal 
places . 

B.A.(sq.  ft.). 

Cross-sectional  area 
in  sq.  ft.  at  top  of 
section.  Three  dec- 
imal places. 

Volume 

— 

-- 

— 

-- 

Cubic  Feet 

Cu.  ft.  volume  compu- 
ted by  method  given  in 
SMALL  BLOCK  item 
"CuC."  (By  Smalian's 
formula  in  present 
computer  program.)  Two 
decimal  places. 

Board  Feet 

— 

-- 

-- 

-- 

Int.  1/4 

Bd.  ft.  volume  ob- 
tained by  International 
1/4"  rule.  One  decimal 
place. 

Scribner 

Bd.  ft.  volume  ob- 
tained by  Scribner 
rule.  One  decimal 
place. 

1/  All  decimal  points  implicit — none  punched  in  IBM  cards. 


2 J Forest  Survey:  PT  is  poletimber-size  tree  from  the  5-inch  D.B.H. 

class  up  to  9.0  inches  D.B.H.  in  Interior  conifers  and  up  to  11.0  other- 
wise. ST  is  sawtimber-size  tree  above  these  limits. 

Timber  Mangt.Res.;  Differs  in  that  minimum  ST  limits  are  8.6-  and 
10.6-inches.  (Computer  Sense  Switch  1 on  for  Survey,  off  for  Timber 
Mangt . ) 

3/  Merchantable  ST  top  taken  as  40%  of  D.B.H.  (outside  bark)  but 
not  less  than  8.0  inches  D.I.B.  in  hardwoods  and  6.0  in  softwoods. 
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Lower  block. — This  irregularly-shaped  block  beneath  the  large  block 
is  for  longhand  entries.  The  computer  program  outputs  similar  quantities. 
Briefly,  the  principal  terms  involved  are  as  follows: 


Term 

Meaning 

d2h 

D.B.H.  x D.B.H,  x Total  height 

F 

Form,  equal  to  FN/FD 

T,  Cu.  ft. 

Cubic-foot  volume,  i.b. , ground  to  tip . 

Ml,  Cu.  ft. 

Cubic-foot  volume,  i.b.,  stump  to 
inch  top. 

4.0- 

Ml,  M/T 

Ratio  Ml,  Cu.  ft./T,  Cu.  ft. 

M2,  Cu.  ft. 

Cubic-foot  vo lume , i.b.,  stump  to 
inch  top. 

6.0- 

M2,  M/T 

Ratio  M2,  Cu.  ft./T,  Cu.  ft. 

M3,  Cu.  ft. 

Cubic -foot  volume,  i.b.,  stump  to 
inch  top. 

8.0- 

M3,  M/T 

Ratio  M3,  Cu.  ft./T,  Cu.  ft. 

Int . , Bd.  ft. 

Board-foot  volume.  International 
inch  rule,  in  16-foot  logs. 

1/4- 

Int. , Bd. /Mr 

Ratio  Int.,  Bd.  ft. /Ml,  Cu.  ft. 

Scrib. , Bd.  ft. 

Board-foot  volume,  Scribner  rule, 
16-foot  logs. 

in 

Scrib. , Bd. /M. 

Ratio  Scrib.,  Bd.  ft. /Ml,  Cu.  ft. 

Cubic-foot  ratios  express 

the  volume  from  stump  to  the  8.0- , 

6.0- , 

-. 0-inch  D.I.B.  tops  as  a 

proportion  of  total  (ground  to  tip) 

vo  lume 

Board-foot/cubic-foot  ratios  express  Scribner  and  International  volumes 
within  the  sawtimber  portion  of  the  bole  as  a proportion  of  cubic-foot 
volume  to  the  4- inch  top.  For  interior  Alaska,  the  cubic-foot  volume  to 
the  8-inch  top  is  figured  from  a 2.0-foot  stump;  cubic-foot  volumes  to  the 
6-  and  4 -inch  tops,  and  board-foot  volume,  are  computed  from  a 1.0-foot 
stump.  For  coastal  Alaska,  all  volume  computations,  except  for  total, 
are  normally  made  from  a stump  height  equal  to  D.B.H.  between  the  fixed 
limits  of  2-  and  4-1/2  feet.  However,  any  stump  height  not  over  4-1/2 
feet  may  be  specified,  since  the  computer  recognizes  it  as  KINDS  01. 

For  both  areas,  the  sawtimber  top  is  taken  as  40  percent  of  D.B.H.  but 
not  under  8.0  inches  D.l.B.  in  hardwoods  and  6.0  inches  D.I.B.  in 
conifers . 

Cross-section  paper. --This  portion  of  the  form  permits  graphic  de- 
termination of  the  cubic-foot  volume.  The  scale  is  10  x 10  to  the  inch. 
Table  1 offers  guides  to  estimate  the  scale  for  the  x-  and  y-axes,  and 
conversion  factors  for  various  combinations  of  maximum  stump  diameters 
and  tree  heights . 
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By  plotting  and  smoothing  old  tree  measurement  data,  one  can  com- 
pute quite  closely  the  diameters  and  cumulative  height  at  any  point  of 
interest.  Entering  such  information  in  the  data  blocks  will  facilitate 
keypunching  and  volume  compilation. 

Computer  Program 

The  FORTRAN  program  was  written  for  the  IBM  1620  having  a 20K  memory 
automatic  division,  indirect  addressing,  and  a 1622  card  read/punch.  The 
debugged  program  and  its  accompanying  block  diagram  (fig.  3)  are  included 
Input  data  cards  are  prepared  as  shown  under  "Tree  Measurement  Form 
Entries  Described"  beginning  on  page  3.  Running  time  per  tree  (averaging 
over  10  input  cards  each)  is  about  one  minute. 

Computer  Output 

Output  by  computer  can  be  typed  and  punched  or  just  punched,  with 
two  typewritten  lines  and  two  cards  per  tree.  Sense  Switch  2 on  the 
console  controls  this  option.  The  first  line  and  output  card  has  all 
the  descriptive  tree  information  and  variables  necessary  for  subsequent 
regression  analysis  by  computer  to  prepare  tree  volume  equations.  The 
second  output  card  for  each  tree  has  supplementary  information  not 
directly  needed  for  later  regression  analysis. 


The  coding  and  meaning  of  the  output  items  for  the  first  typed  line 
or  punched  card  follow: 


Code 

NOTRE 

Ell 

DBH 

TLGTH 

F 

HTLOG 

VSCR 

VINT 

VM3 

VM2 

VM1 

VT 


Meaning 

Tree  number 

Sub-region,  unit,  and  species. 

D.B.H.  or  D.O.B.  above  swell. 

Total  cumulative  length,  ground  to  tip. 
Form 

Number  logs  to  closest  16-foot  log. 

Scribner  board-foot  volume. 

International  1/4-inch  board-foot  volume. 

Cubic-foot  volume  inside  bark  to  an  8.0- 
inch  D.I.B.  top. 

Cubic-foot  volume  inside  bark  to  a 6.0- 
inch  D.I.B.  top. 

Cubic-foot  volume  inside  bark  to  a 4.0- 
inch  D.I.B.  top. 

Cubic-foot  volume,  inside  bark,  ground 
to  tip. 
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The  coding  and  meaning  of  the  output  items  for  the  second  typed  line 
or  punched  card  follow: 

Code Meaning 

NOTRE  As  above 

Ell  As  above 

EI2  Age  at  stump  and  local  forest  type 

EI22  Elevation  in  hundreds  of  feet 

EI3  Form  class  code,  cubic-foot  volume 

determination  code,  and  board-foot 
volume  determination  code. 

RVSCR  Ratio  of  VSCR  to  VMl 

RVINT  Ratio  of  VINT  to  VMl 


RVM3 

Ratio 

of 

VM3 

to 

VT 

RVM2 

Ratio 

of 

VM2 

to 

VT 

RVM1 

Ratio 

of 

VMl 

to 

VT 
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Figure  1. — Tree  measurement  form  developed  by  the  Northern  Forest  Experiment  Station  showing  sample  entries 


Table  1. --Guide  to  estimate  intervals  for  X-  and  Y-axes  and  conversion 

factors  for  graphic  determination  of  cubic  volume  for  specific 
maximum  stump  diameters  and  tree  heights 


Maximum  d.i.b.  :statis-:  Maximum  total  height  of  tree  (feet) 


(inches) 

tic_w 

50 

100 

150 

200 

y 250 

Y 

4 

5.6 

X 

5 

Cf 

0.001091 

Y 

10 

10 

8.9 

X 

5 

10 

Cf 

.002727 

0.005454 

Y 

25 

25 

14.0 

X 

5 

10 

Cf 

.006818 

.013635 

Y 

50 

50 

20.0 

X 

10 

15 

Cf 

.027270 

0.040905 

Y 

100 

100 

100 

28.2 

X 

10 

15 

20 

Cf 

.05454 

.08181 

0.10908 

Y 

200 

200 

40.0 

X 

15 

20 

Cf 

.16362 

.21816 

Y 

300 

300 

48.9 

X 

20 

25 

Cf 

.32724 

0.40905 

Y 

400 

i/56.5 

X 

25 

Cf 

.54540 

i/63.2 

Y 

500 

X 

25 

Cf 

.68175 

Y 

1,000 

i/89.4 

X 

25 

Cf 

1.36350 

1/  Y,  X — intervals  for  the  one-inch  points  on  Y-  and  X-axes. 

Cf  — conversion  factor  to  multiply  times  number  of  0.1  inch 
squares  beneath  curve,  giving  cubic  foot  volume. 

2_/  Hand  compute  or  use  larger  graph. 
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Figure  2. --Li  sting  of  FORTRAN  program. 

C THIS  1620  FORTRAN  PROGRAM  READS  INFO  PUNCHED  FROM  TREE  MEAS.  FORMS. 

C 

C FOR  INTERIOR  ALASKA  SPP.  , COMPUTES  CU.  FT.  INSIDE  BARK  BY  SMALIANS  FOR 
C GROUND  TO  TIP,1  FT  STUMP  TO  4 . 0 AND  6.0  INCH  TOPS,  AND  2 FT  STUMP  TO  8 
C .0  IN.  TOP.  ALSO,  SCRIBNER  AND  INT.  1/4  INCH  RULE  VOL.  IN  16  FT  LOGS— 

C FROM  1 FT  STUMP  TO  40  PERCENT  OF  DBH  (BUT  NOT  UNDER  8.0  IN.  TOP  IN 
C HDWDS.  OR  6.0  IN.  TOP  IN  SFTWDS) 

C 

C FOR  COASTAL  ALASKA  SPP.  , COMPUTES  CU  FT  BY  SMALIANS  FOR  MERCH.  STUMP  TO 
C 8.0, 6. 0,4.0,  AND  GROUND  TO  TIP.  ALSO,  INT.  1/4  INCH  AND  SCRIBNER  FROM 
C MERCH.  STUMP  TO  TOP  DESCRIBED  ABOVE. 

C 

C OTHER  OUTPUT  IS  TREE  NO.  , AREA,  SPP.  , DBH,  T.  HT,  FORM  , 3 CU  FT  VOL.  RATIOS, 

C BD  FT/CU  FT  RATIOS, AND  NO.  OF  LOGS. 

C OUTPUT— 2 PRINTED  LINES  AND  2 CARDS  PER  TREE. 

C 

C PAUL  M HAACK,  NOR. FOR.  EXP.  STA.  JUNEAU, ALASKA, REV.  12/31/62. 

C 

C OPERATOR- --CLEAR  MEMORY  AND  TABS, AFTER  ENUF  PRINTED  ANSWERS  ARE 
C VERIFIED,  HIT  READER  STOP,  AND  WHEN  TYPING  STOPS,  TURN  OFF  SENSE 
C SWITCH  2,  HIT  READ  START, AND  PROGRAM  WILL  CONTINUE  PUNCHING  ONLY. 

C SWITCH  1 OFF  FOR  TIMBER  MANGT.  ST  MINIMUM, ON  FOR  FOREST  SURVEY. 
IF(SENSE  SWITCH  9)67,67 
67  IDTRE=0 
TVM  2=0.0 
TVM  3=0.0 
RVM1=0 . 0 
RVM  2=0.0 
RVM  3=0.0 
RVINT=0 . 0 
RVSCR=Q . 0 
VT=0 . 0 
VM  1=0.0 
VM2=0 . 0 
VM  3=0.0 
VINT=0 . 0 
VSCR=0 . 0 
CLGTH=0 . 0 
TLGTH=0 . 0 
HTLOG=0 . 0 
INO  S=0 
ADTO3=0 . 0 
2 INOS=INOS+l 

READ100,  NOTRE,  Ell  ,EI2,EI3,  DBH,  HT,FN,FD,  NOS, KINDS,  SLGTH,  DBT,  DIB 
100  FORMAT(l  3 , F5 . 0 , F5 . 0 , F8 . 0 , F3 . 1 , F5 . 2 , F3 . 1 , F3 . 1 , 12 , 12  , F4 . 2 , F3 . 1 , F3 . 1) 
IF(IDTRE)69 ,70,69 

(continued) 
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70  IDTRE=NOTRE 
69  IF(IDTRE- NOTRE)  43, 23,43 
23  IF(INOS-NOS)  802,230,802 
230  GO  TO  (3, 5, 5, 5, 5, 7, 7, 7, 7),  KINDS 
3 VT=VT  +(DIB*DIB*0.005454*SLGTH) 

STUMP=SLGTH 
BIGDI=DIB 
GO  TO  2 

5 AVEBA=(0.005454*(  (BIGDI*BIGDI)+(DIB*DIB)  ) ) /2 . 0 
SECV=AVEBA*  SLGTH 

VM1=VM1+SECV 
BI  GDI=DIB 

CLGTH=CLGTH  +SLGTH 
IF(KINDS-2)  801,6,8 

6 ADT03=VM1 
GO  TO  2 

8 IF(30000 . 0 -Ell)  56,56,57 

56  IF(SENSE  SWITCHl)352 , 58 

58  IF(DBH-8. 6)320, 354,354 

57  IF(SENSE  SWITCH1)353 , 59 

59  IF(DBH-10. 6)320, 10, 10 

352  IF(DBH-9.0)  320,354,354 
354  T SP=EI1*100 . 0 

IE(ISP-6900)10, 10,57 

353  IF(DBH-11 .0)320, 10, 10 
10  IF(DIB-8 .0)51 , 50,110 

50  TVM3=VM1 
GO  TO  110 

51  IF  (DIB- 6 . 0)54 , 54 , 1 10 
54  TVM2-VM1 

110  IF(CLGTH-16 . 3)2 , 12 , 12 
320  TLGTH=TLGTH  +CLGTH 
CLGTH=0 . 0 
GO  TO  2 

12  VSCR=VSCR+(  (0 . 79*DIB*DIB)-(2 . 0*DIB)-4 . 0) 
VTNT=VINT+(  (0 . 796*DIB*DIB)-(1  . 37  5*DIB)-1  .23) 
HTLOG=HTLOG  +1.0 
TLGTH  =TLGTH  +CLGTH 
CLGTH=0 . 0 
GO  TO  2 

7 IF(CLGTH-8.15)  28,30,30 

30  VSCR=VSCR+(  (0 . 79*BIGDI*BIGDI)  - (2.0*BIGDI)  -4.0) 
VINT=VINT+(  (0 . 79  6*BIGDI*BIGDI)-(1 . 375*BiGDI)-l . 23) 
HTLOG=HTLOG  +1.0 
28  TLGTH=TLGTH  +SLGTH  +CLGTH 
CLGTH  =0.0 

AVEBA=(0 . 0054 54* ( (BIGDI*BIGDI)+(DIB*DIB)  ) )/2.0 
SECV=AVEBA*  SLGTH 

GO  TO  (800,800,800,800,800,32,34,36,38),  KINDS 


(continued) 
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32  VM 1 =VM  1 +SECV 

VM3=VM1 
GO  TO  360 
34  VM1-VM1  + SECV 
VM2=VM1 
GO  TO  360 
36  VM1-VM1  + SECV 
360  BIGDI=DIB 
GO  TO  2 

38  VT=VT+(  (BIGDI*BIGDI*0.005454*SLGTH)/2.0)+VMl+0.05 
TLGTH=TLGTH+STUM  P 
IF(HT-TLGTH)40, 20,41 

40  IF(HT-(TLGTH-1 .0)  ) 803,20,20 

41  IF(HT-(TLGTH+1 .0)  ) 20,20,803 
20  VM1=VM1  +0.05 

IF(VM1  -0.05)358,358,359 
358  VM1=0 . 0 
RVM 1=  0 . 0 


GO  TO  36  2 

359  RVMl=(VMl/VT)+0 . 0005 

362  IF(TVM2)  364,364,363 

363  VM2=TVM2  +0.05 
GO  TO  367 

364  IF(VM2)  366,366,365 

365  VM2-VM2  +0.05 

367  RVM 2 =(VM2/VT)+0 .0005 

366  IF(TVM3)  370,370,369 

369  VM3=TVM3  +0.05  -ADT03 
GO  TO  373 

370  IF(VM3)  372,372,371 

371  VM3  «VM3  +0.05  -ADT03 
373  RVM3=(VM3/VT)+0 . 0005 

372  IF  (VINT)  4 2 ,42,355 
355  VINT=VINT  +0.5 

RVINT=( VINT/VM 1 )+0 . 0 0 0 5 
VSCR  =VSCR  +0.5 
RVSCR=(VSCR/VMl)+0 . 0005 
42  F=FN/FD 

IE3-EI3/1 000000. 

EI22=IE3 

DI22-EI2 2*1 000000. 

EI3=EI3-DI22 

IF(SENSE  SWITCH  2)  208,209 

208  PRINT1 , NOTRE,  Ell , DBH,  TLGTH,  F,  HTLOG,  VSCR,  VINT,  VMS,  VM2,  VM1 , VT 
1 FORMAT*  14,  F7.0,F5.1  ,F7  .2  , F5  .3  ,F4  .0  ,F6  .0  ,F6  .0  ,F6  .1  ,F6.1  ,F6.1  ,F6.1) 

PRI  NT  2 0 1 , NO  TRE , El  1 , El  2 , El  2 2 , El 3 , RVS CR , RVI  NT , RVM  3 , RVM  2 , RVM  1 
201  FORMAT(  14  , F7 . 0 , F7 . 0 , F4 . 0 , F8 . 0 , F6 . 3 , F6 . 3 , F6 . 3 , F6 . 3 , F6 . 3//) 

209  PUNCH1 , NOTRE,  Ell , DBH,  TLGTH,  F,  HTLOG,  VSCR,  VINT,  VM3,  VM2,VM1 , VT 
PUNCH  201 , NOTRE,  Ell,  EI2,  EI22,  EI3,  RVSCR,  RVINT,  RVM  3,  RVM  2,  RVM1 


IF(SENSE  SWITCH  9)  206,67 
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206  PRINT  207 

207  FORMAT(48HYOU  HIT  READ  START  TO  GET  LAST  CARD,  PROGRAM  DONE//) 
STOP 

43  PRINT  202 

202  FORMAT(4 1 HCHECK  TREE  NO CORRECT, REPLACE,  PUSH  START//) 

PAUSE 

GO  TO  23 

800  STOP 

801  STOP 

802  PRINT  203 

203  FORMAT(44 HCHECK  SECTION  NO.  , CORRECT, REPLACE,  PUSH  START//) 
PAUSE 

GO  TO  230 

803  PRINT  204 

204  FORMAT(46HCUMUL.  HT.  OFF  OVER  FT.  FROM  HT,  MAYBE  RERUN  TREE//) 
PAUSE 

GO  TO  20 
END 


(THREE  OUTPUT  EXAMPLES) 

+1+11448. +21 .6+99 .70+. 800+3.  +368.  +405.  +68.6  +72.7  +75.0  +80.5 
+1+11448. +18048.  +5.  +30103.  +4.907+5.402  +.852  +.904  +.932 

+35+33144.  +8.9  +50.00+.820  +.  +.  +.  +.0  +8.7  +10.1  +11.0 

+35+33144.  +7041. +10.  +30100.  +.000  +.000  +.000  +.792  +.922 

+ 36+33144.  +9 .0  +50. 10+. 811  +2.  +36.  +50.  +.0  +8.8  +10.2  +11.1 
+36+33144.  +7041. +10.  +301 01 . +3 . 540+4 . 9 38  +.000  +.794  +.923 

YOU  HIT  READ  START  TO  GET  LAST  CARD,  PROGRAM  DONE 


STOP 
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Diagram  of  FORTRAN  Program. 
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